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ABSTRACT 

Electromagnetic waves impinging upon a plasma at frequencies larger than 
the plasma frequency, suffer weak scattering. The scattering arises from fhe 
existence of electron density fluctuations. The so-called incoherent scattering 
tiieory basically deals with fluctuations of random thermal origin; however, for 
practical purposes, it must also take account of those fluctuations caused by 
streaming photo-electrons. As is well known, in any scattering experiment, the 
received signal corresponds to a particular spatial Fourier component of the 
fluctuations, the wave vector of which is a function of the wavelength of the radio- 
wave. Wavelengths short with respect to the Debye length of the medium relate 
to fluctuations due to non-interacting Maxwellian electrons, while larger wave- 
lengths relate to fluctuations due to collective Coulomb interactions. In the latter 
case, the scattered signal exhibits a spectral distribution which is characteristic 
of the main properties of the electron and ion gases and, therefore, provides a 
powerful diagnosis of the state of the plasma, in our case, the ionosphere. 


*NAS/NRC Senior Research Associate, On leave from the Centre National d’Etudes des 
Telcommunications, Department RSR (France) 
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INTRODUCTION 


The steps leading to the theory of elecjtrom^netic waves incoherently 
scattered by the ionosphere, or more generally by a plasma, for operating 
frequencies much larger than the plasma frequency, can be tracked back to the 
last century. It was not until 1958, however, vhen Gfordon (1958) became aware 
of the fact that such a scattering was within reach of modern radars, that a large 
and very successful effort was conducted in order to establish a complete theory. 

The first keystone mi^t be thou^t to be the theory of the scattering of 
light by the molecules of the atmosphere, developed by Lord Rayleigh (1871). 

The next and major step is the work of Thomson (1906) who showed that a free 
electron acting as a dipole is capable of scattering electromj^etic waves at all 
frequencies. Therefore electromagnetic waves ought to suffer a small scattering 
from a plasma in a maimer similar to the scatterir^ of light by molecules. This 
last point was further emphasized by Fabry (1928) who indicated that the random 
motion of the electrons should induce a doppler broadening of the back-scattered 
signals. The importance of the random motion of the electrons on the scattering 
has later led to calling the phenomenon 'Tncoherent scatter" rather than "Thomson 
scatter". A large effort followed to develop the theory of random motions in 
plasmas and credit should be given to cite but a few, to Landau (1946), Bohn and 
Gross (1949a, 1949b), and Pines and Bohn (1952). 

Interest in ionospheric research developed rapidly after pioneering experi- 
ments of Appleton and Barnett (1925), Breit and Tuve (1925), and Taylor and 
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Hulburt (1926) -wMch definitely located the ionospheric layer. By 1931 Chapman 
gave a theory of the ionizing effect of monochromatic radiation in an atmosphere 
on a rotating earth. The available techniques however put a limit to the study 
of the ionosphere to essentially the bottom side electron density profiles, leaving 
out, for example, experimental studies of the energetics of the ionosphere. 

The theory of scattering of electromagnetic waves by the earth environment 
was also in progress. Booker and Gordon (1950) considered radio scattering 
from the troposphere and Booker (1955) scattering by nonisotropic irregularities 
in the aurora. 

Gordon (1958) finally took the decisive step vhen he realized that the in- 
coherent scatter from free electrons of the ionosphere would amount, provided 
powerfiil radars be used, to a sufficiently large fraction of the sky noise, even 
taking into account the fact that the doppler broadening due to the thermal motion 
of the electrons was expected to be quite large. The total cross section was 
computed to be of the order of .25 cm^ for a i" antenna beamwidth and a 67^xs 
pulse length with a mean electron density of 10^ cm ^ (Farley 1970). Gordon 
pointed out that this experiment would provide electron density and temperature 
data up to several earth radii. 

The suggestion made by Gordon met very favorable circumstances. Indeed, 
it coincided with the advent of the space age, that is, as far as the ionosphere is 
concerned, with direct probing of the medium and consequently an increasing 
interest in the understanding of the processes governing it. 
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The first incoherent scatter experiment was performed the same year by 
Bowles (1958) using an already existing transmitter together with an array of 
dipoles. The return power appeared to be ^proximately of the expected order 
of magnitude but the bandwidth was observed to be much narrower thgr> predicted 
by Gordon. 

While several groups were in the process of developing incoherent scatter 
facilities for the purpose of ionospheric studies, several authors undertook the 
studies of the electron Ouctuations including the effect of the ions which rapidly 
appeared to be responsible for the narrow power spectrum observed by Bowles. 
The case of independent electrons would appear as a limiting case, and in most 
instances the small degree of organization of the medium, due in particular to 
Coulomb interactions, ought to dominate the characteristics of the scattering 
process. The theory of the scattering had become much more involved than 
thought in the first place, in turn a much more detailed description of the medium 
was to be obtained from the incoherent scatter spectrum. 

The scope of this paper is to outline the different steps of the building of 
the theory. The geophysical results as well as the technical description of the 
existing facilities, which have been extensively reviewed by Evans (1969, 1974) 
and Earley (1970b), are not covered here. 

I. Thomson Scattering of an Electromagnetic Wave by a Free Electron 

When a field E^- is applied to a free electron at rest, an 

oscillatory motion is Imparted to such as: 



1.1 


V = 






m 


where ojq is the angular frequency, t the time, q the charge electron, the 
electron mass and m, ^ the time rate of change of momentum of the electron. 
In turn the acceleration of the electron gives a rise to a radiation field E, . 


E = 


1 R X (R X V) e 


1.2 




where R is the radius vector joining the position of the electron to the field 
point, c the speed of light in vacuum, dielectric constant of the vacuum. 

Use of I.l yields: 


E. = 


„2 Rx (RxE.) -(i<iJ R/c) 

^ e 
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therefore the magnitude of the field is 


3 ^ ^ 


E„ = 


4776„C^m ^ 

O e 


1.4 


where ^ is the angle between Ep andR. The same treatment applies to the 
ions but because of the mass difference the radiation field is negligible compared 
to that of the electrons. Equation 1.4 can be written as: 


E = En ^ 

s JJ 0 




sin ^ 


1.5 
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where ) is the classical electron 

radius: r =2.82 x 

e 

More generally R is a function of time, and the frequency of the scattered field 
is doppler shifted with respect to the firequency of the incident field Ej . 

The mean radiated power per unit area is 


P 




s in^ ^ 


1.6 


where is the m^netic permeability in vacuum. Equation 1.6 can be written 
as 



R2 


1 

2 
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where = r^ sin^ ^ Is the differential cross section. 

II. Thomson Scattering by hidependent Randomly Distributed Free Electrons 
The problem is now extended to a set of independent randomly distributed 

I 

free electrons. 

The effect of multiple scattering can be neglected in practical cases because 
of the smallness of the differential cross section. The working frequency is 
supposed to be much larger than the frequency of free oscillations (plasma 
frequency). It follows from the above that the transmitted wave is not signifi- 
cantly altered: The Born approximation is valid. It is farther assumed that 
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the electrons of interest are confined in a small region compared to a sphere 
of radius R, where R is the mean distance between the field point and the scatter' 
ing electrons. 

The total field at the field point can be e 3 q>ressed as 

II.2 


E sin^E^e ^ e 


'^i 


n.2 


where j indicates a summation over all the electrons and 0. is the phase 
associated with the location of each of the electrons relative to both the trans- 
mitter and the receiver sites. However, because of the random nature of the 
electron distribution in space, the instantaneous field at the field point is of no 
interest. The average radiated power per unit area is the physically interesting 
quantity. 
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where Rg desigjiates the real part and the brackets indicate the ensemble 
average. The first term in the square brackets is just the number of 
electrons in the volume of interest. The second term represents the mean 
phase relationships of the electrons taken two by two; this term has to be zero 
since the electrons are supposed to be independent. 

The radiated power is therefore simply the sum of the radiated power of 
each electron taken independently 

<p>=2^Pi 

j 

and the total differential cross section is 

cr. * = N Cr^ 

tot e e 

This is the result obtained by Gordon (1958). As for the doppler broadening of 
the scattered field, it simply corresponds to the velocity distribution of the 
electrons. In the case of backscatter the doppler shift associated with a velocity 
V is (o) -Wq) = -2w(jv/c. Therefore to a Maxwellian distribution of velocities 


f (v) a 


- m V 


2KT 


(where K is the Boltzmann constant and Te the electron temperature) there 
corresponds a power spectrum 


S(oi) ae 


-m^/2KT„ , 1 

®(c(o;-a>o)/2<yo)^ 


n.7 
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The first experiment of Bowles (1958), however, showed that in general the 
neglect of the term 



is not valid. This term represents the degree of organization of the medium 
which arises from particle interactions. 

The object of the rest of this paper is to investigate the degree of organiza- 
tion brought into the medium under various circumstances. It will be shown that 
in general the total differential cross section obtained for the case of independent 
electrons is an upper limit. 


in. Formulation of the Scattering Problem 

It is Srst convenient to express the scattered field in the following way 


E 


S 



in.i 


where is the wave vector of the transmitted electric field, k 2 the wave vector 
of the scattered field in the direction of observation and f j the radius vector 
between the j electron and the center of the scattering volume (the incident 
field at the center of the scattering volume is just E<, e ). 

De fining k as - k^ leads to 


E = — s 
= R 


V » • V. ) 

ic ^^0 / ^ g 


m.2 
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for point charges the electron distribution can be expressed as; 



j 


the spatial fourier series of which is 

N (r, t) N (k, t)|^ 

k 

such that 


in.3 


in.4 




where V is the scattering volume. Upon substituting Eq. III.5 into IH,2 yields 


Es(t) = — (k, t)|^ 

R 


which means that the scattering only arises from electron density fluctuations 
and more precisely from the k fourier component of these fluctuations. Stated 
in other words the scattering only depends upon stratified fluctuations of the 
density perpendicular to the k direction (Fig. 1). 

Equation ni,6 also shows that the scattering properties to be investigated 
depend upon the observing scale 1 in fact it will turn out that the degree of 
organization of the fluxtuations of the medium responsible for the scattering 
depend upon the relative magnitudes of 1/k and the Debye length of the plasma 
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\.= 



witix No the electron number density. The averaged scattered power is again 
given by 


m.7 

= <^,(k, t)l N;(k. t)| \ 

2 r /Xq p2 


By virtue of the Wiener-Khintchine theorem the power spectrum is just the 
Fourier transform of the autocorrelation function; 



in.8 


Thus the power spectrum is then given by 

PsH + ^) =4]/^ <1E>0 + ^)l' > 

^ m.9 
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In the following the only quantity of interest is consequently 


<(M. 

The associated total differential cross section 

+ “) = Te ^ ^ <^Ng(k, in.lO 

The determination of (k, cu) , the spectrum of the random electron 
density fluctuations, requires in general the solution of Maxwell equations and 
of the kinetic equations for the perturbed plasma. 

The various parameters which affect the spectrum are the Coulomb inter- 
actions, the magnetic field, the collisions with neutrals, the electron and ion 
temperatures, the drifts of the ions and of the electrons, the ion composition 
and the presence of a non maxwellian tail in the electron distribution. The main 
difficulty is to take account of these by using the proper kinetic equations. 

For each situation a different dispersion relation for the longitudinal plasma 
waves is associated with the scattering problem. The wave modes and the 
different damping processes involved will be of help for the interpretation of 
the shape of the spectrums. 

The solution of the different sets of equations also requires a statistical 
treatment because of the random nature of the initial electron distribution in 
space. The analysis performed by a number of authors differs mostly in the 
way the statistical treatment is carried out; Fejer (1960, 1961), Salpeter (1960), 
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Renau (1960), Hagfors (1961), considered the distribution probability of the 
electrons along lines similar to the one followed Pines and Bohm (1952) while 
Dou^erthy and Farley (1960), (1963) Farley et al. (1961) made use of the 
Nyquist theorem. The identity of the results obtained by the different methods 
indicates their consistency. For the sake of mathematical simplicity the 
approach rn ptdng use of the Nyquist theorem is to be outlined in this paper. 

It must be noted (Farley (1970)) that the use of kinetic equations can be 
avoided for two classes of problems; 

Firstly, if just the total scattered power is needed, use can be made, in 
some cases, of the Debye-Huckel theory (Kahn 1959; Salpeter 1960; Renau 1960). 
In this connection, see the proceeding section in which the total power scattered 
by independent electrons was derived without any assumption about the electron 
velocily distribution. 

Secondly, fluid equations can be used instead of kinetic equations in the 
case of collision-dominated plasmas (Tanenbaum 1968, Seasholt and Tanenbaum 
1969), 

m. summary, the scattering problem involves one particular Fourier 
component of the electron density fluctuations depending upon the relative 
locations of the transmitter, the scattering volume and the field point. The 
determination of the power spectrum of this Fourier component implies the 
use of Maxwell’s and kinetic equations in the (k, w ) space. The shape of the 
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spectrum can in general be simply interpreted in terms of heavily damped 
longitudinal plasma waves. 

rv. The Nyquist Theorem Approach 

The procedure followed by Dougherty and Farley (1960) for obtaining 
<|N, (k, >i:> is outlined in this section. 

The Nyquist theorem in its most simple form states that in an isolated 
resistance R at a temperature T there flows a random current I such that 

<|i(.)|’>do.=i IV.1 

^ 7T K 

The theorem can be extended to a linear system in thermodynamic equi- 
librium at a temperature T whose response I(^) to a generalized force V(^ )is 

I(co) = Y(co) V(o;) IV.2 

where Y is the general admittance function. In the absence of the driving 
force the response of the system is such that; 


<f|l(£o)r'>doj = iG(co) KTdw IV.3 

^ ' 77 

where 


G(co) = (Y(o>)) 

and o) varies between-® and +®. 
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When applied to a plasma the theorem is often called the fluctuation- 
dissipation theorem, whereby random thermal fluctuations in the electron 
density can be related to a dissipation process such as Landau damping or 
collisions (Farley, 1970a). 

The way of applying the theorem which is relevant to this problem is to 
calculate the displacement of the electrons in response to a fictitious force f 
applied to them only. However the motion of the electrons and of the ions are 
coupled through the space- charge electric field E. Use of Maxwell's equations 
allows for the determination of the generalized admittance tensor Y ' in terms 
of the admittance tensors of the electrons and ions Yg and Yj • 

The generalized admittance, in turn, leads to the spectrum of the density 
fluctuations. The determination of the admittances Y. and Y calls for the 
kinetic equations of the electrons and ions. This will be the object of the next 
sections for various conditions. Assuming charge neutrality and singly charged 
ions the mean flow vector resultii^ from f are; 


N .U = _ ^ 

■ (F + qE) 

V ® 


IV.5 
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where N^o is the unperturbed electron number in the volume V, and and 
the electron and ion drift vectors. Introducii]^ the current 

J=^q(U,-U,) IV.6 

and eliminating B in the MaxweD’s equations written in the (Ic, co) space yields; 


J = - i€.E + — X (kx E) 


IV.7 


or 


J = r • E 


rv.8 


r is a diagonal tensor if a coordinate system is chosen in such a way one 
axis is parallel to k. Eliminating E, J and Uj between equations IV.4, IV.5, 
IV.6 and rv.8 yields 


N,oU, = VY'-F 


IV.9 


with 



Use of the continuity equation of the perturbed electron density n^ yields 



= 0 


IV.11 
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which can be linearized in (k, o)) space: 


N ^ ^ 

VN^ (k, w)l = - k • U^(k, 0)) 
® y CO 


rv.i2 


therefore 




k2 




<lNeoUL(k. 


IV .13 


and by virtue of the Nyquist theorem 

V= < iN.Ct ")i:>= 4 -)1^ = ktvR.(y;,) 

^ ® TTO/ 


finally 




IV .15 


V. The Case of Collisionless Plasma in Thermal Equilibrium without 
Magnetic Field 

The computation of the admittances Yi and Ye for the ions and the electrons, 
which enter in (IV .15), is performed using the appropriate Boltzmann equation 
and a perturbation method. The electron velocity distribution can be written 

as; 

f (v) = 
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where is a Maxwellian distribution and is a perturbation. 

r -|3/2 

r(^) = nJ-!!L_ e--”’ /2KT V.2 

® ® L27tKTJ 

\\iiere m is the mass and no the number density of the particle of interest, 

— ♦ 

The Boltzmann equation, to the first order in perturbations in (k, co) space is: 

i(o>-k'v)f, + i.F* — ^ =0 V.3 

' 

The next step consists in evaluating n<j U 

• U = J vfj d^v = — F * / — 2 d^v V.4 

^ dw CO -k'v 

OqU = f • F V.5 

Assuming again that k is parallel to one of the coordinate axes it can easily 
be shown (Dougherty and Farley 1960) that Y is diagonal. Therefore all the 
tensors entering in Y' are diagonal and only the kk element of each needs to be 
considered. 

Before givir^ a final expression for Y it is interesting to recall the physical 
meaning of the singularity in V.4 for co = k • v . 

hi terms of waves, the singularity corresponds to particles movu^ at the 
phase velocity of the wave. Consequently they tend to see a constant accelerating 
or decelerating electric field, and interact strongly with the wave. Ultimately 
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after a net gain or loss of energy the particles oscillate in a reference frame 
moving at the phase velocity. The linear approach therefore breaks down since 
it predicts a steady exchange of energy between the particle and the wave. This 
effect has been considered by Landau (1946) and is known as "Landau damping". 
It can also be shown (see Bohn and Gross, 1949, for example) that in addition 
to the solutions {a > , k) of the dispersion relation for a collisionless gas of one 
type of particles, wave energy can be associated with each value of w for which 
there exist particles whose velocities satisfy the relationship k • v = co . 

It can also be inferred that the power spectrum of the fluctuations corres- 
ponding to these waves decreases when k decreases, while the opposite is true 
for the waves predicted from the dispersion relation and for which the contri- 
butions of all particles add coherently. 

The mathematical difficulty implied by the singularity was overcome by 
Landau (1946) and the result is (Dougherty and Farley (1960): 



where 




1/2 


or 


‘kk 


= — 

mct) 


V.7 
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where 


W(0) = 1 - 



dp + in 



Therefore substituting (0^) and Y^y. (0^) into eq* IV .10 yields 


^kk 


in^o; [a^W(<9, ) + 1] 

^ 

k^KT [a.2(W(5i) + W(^^)) + 1] 


with 


V.8 


V.9 


k^KTe^ k2\2 

a relates the scale of the observation, determined by the geometry of the ex- 
periment and by the wavelength, to the characteristic scale of the plasma, 
determined by the plasma number density and the plasma temperature. There- 
fore the power spectrum of the density fluctuations is: 




-2- s 

7TO)V ® 


iW(0J [a2W(0.) + 1] 
Ja2(W(^,) + W(<9J) + 1]^ 


V.10 


or 




V/2ifV[2KT]^^2 


V.ll 






-02 


la2w(^^) + l\ + mpe " | a2W((9^)| 


a2[W(^.) + W(0J] + 1 
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Noting that the \W(0) and 3^W(0) are of the order of 1 or smaller, two limiting 
Cases of practical interest can be considered, according to the relative magni“ 
tudes of the observing and plasma scales, 
a) a < < 1 

Equation V.ll becomes 




m 


-il/2 


-<u^nig/k^2KT 




L2KTJ 


V.12 


which is just the result of Gordon (1958): Therefore the spectrum simply 
corresponds to independent free electrons. The physical interpretation of this 
situation is that no collective oscillations can be sustained on a scale smaller 
than the Debye length, 
b) a > > 1 

It is first interesting to note that equatii^ to zero the denominator of eqn. V.11 
yields the dispersion relation for plasma oscUlations (Landau 1946): 

a2(W(^.) + W(5J) + 1 = 0 V.13 

the solution of Eqn. V.13 gives two modes: the first is the pseudo acoustic 
mode which can be viewed as sound waves of a gaz consisting of particles having 
the mass of the ions and the sum of the energies of both the electrons and the 
ions. The second corresponds to electrostatic waves at a frequency close to 
the plasma frequency. 
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Because of the mass difference between ions and electrons the numerator 
of eqn. V ,n, is dominated by the second term when ^ 1 and by the first term 
when 1. The corresponding simplified expressions are termed respectively 

ionic and electronic spectrum. 
b.l) The ionic spectrum 1» 6^ <<1) 

Eqn. V .n. reduces to: 




n„ 





-a? 
e * 


![ w ( 0 .) + 1]|2 


V.14 


This shows that the spectrum corresponds to an ion velocity distribution modi- 
fied by the term in the denominator. Figure 2 (after Farley 1970) is a diggr ^tnn 
of the normalized differential cross section for backscatter per unit volume 
cr (0. )7T^''2/{nn, r^) where a- {6.^) is defined by the relation 


+ O)) = 1 ^ VE2<t(<9,) de, 


V.15 


Also shown on the diagram is the normalized ion velocity distribution e ^ . 

Such a spectrum was first observed by Bowles (1958). From the width 
of the spectrum be concluded that the ions are playing a major role. The influ- 
ence of the denominator affects however significantly the spectrum which is 
depressed in the center and enhanced on the sides because of the presence of 
heavily damped pseudo acoustic waves. 
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The damping of the waves reflects the fact that the wave phase velocity 
lies in a range of the ion velocity distribution such that many particles interact 

with the wave and therefore cause daoiplng* 

It might be worth noticing feat at the time the first incoherent scatter 

measurements were made, the ions were generally considered as a uniform 
charged background. 

The measurements therefore provided an immediate experimental support 
to new theories of plasma oscillations taking into account the discreet nature 
of the ions such as the one of Kahn (1959). Incoherent scatter observations 
therefore played a significant part towards the understanding of plasma 

, oscillations. 

RecalUag that the electrons are responsible for the scatter It is interesting 
to give an Interpretation of the reason why the observed characteristics match 
closely these of the ions# 

If the ions were by themselves thermal oscillations on a scale much larger 
than the Debye length (that is a > > 1), would be highly restricted because 
of the opposing effect of the space charge electric field which would arise from 
those fluctuations. The conditions are however quite different because of the 
presence of the electrons. The electrons which have a large mobility compared 
to that of ions respond very effectively, to the space charge which tends to develop 
and this has the following effects; (i) the electron fluctuations nearly match the 
ion fluctuations, and therefore the spectrum of the electron fluctuations pictures 
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the ion fluctuations, (ii) the space charge electric field is considerably reduced 
with respect to the one obtained from the same fluctuations in the absence of the 
electrons. The residual electric field is comparable to the ambipolar (polar- 
ization) electric field and leads to pseudo neutral having the mass of the ions 
and the sum of the energies of the electrons and the ions, (iii) the residual 
electric field, however small, is sufficient to distort significantly the spectrum 
of the ion fluctuations with respect to the one obtained for non interacting ions. 
b,2) The Electron Spectrum 6 ^ > 1 , 6 ^ >> 1 
Equation V.l reduces to 




k7^1/2v 



+ 1 |^ 

a2W(6'J + l|2 


V.16 


This is generally a small quantity except for values of 6 ^ for which the de- 
nominator is a minimum, that is for the solutions of the dispersion relation 
corresponding to electrostatic waves at about the plasma frequency. 

In the corresponding ranges of 



equation V.II further reduces to 



V.17 
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The resonant pulsation is given by {Bohm and Gross 1949 for example): 


O? + 

T p ^ 



V.18 


where Wp corresponds to the plasma frequency, or in normalized units: 

02 3 ?!..+ ,?. V.19 

er 2 

the corresponding spectrum consists in two very narrow lines called the plasma 
lines; their width is, according to the initial assumptions, determined only by 
the Landau damping: 




1/2 4 — 0,^/2 

TT^' a’ e 


2 


V.20 


The Landau damping being small, because of the small number of electrons 
moving at the phase velocity of the waves, other processes such as collisions, 
not taken into account here might play a significant role. This will be investi- 
gated in other sections. 

c) Figure 3 from Hagfors (1961) shows complete normalized spectra for various 
values of a . 

For a/ i/T = 300 the spectrum has a width corresponding to the ion velocity 
distribution. The plasma line, missing in this diagram, would be found for 
0 - 300; however the corresponding power would be negligible. Therefore the 
spectrum is essentially an ion spectrum, for a/v^ = .1 the spectrum corresponds 
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to independent free electrons. In particular the plasma line which should be 
found for 6 ~ is completely damped because the corresponding phase 

velocity is close to the thermal velocity. 

Intermediate cases such as a / 3 exhibit, with different amplitudes, 
both the electronic and ionic velocity distributions characteristics; the pla s ma 
line is also present but contributes to a negligibly small part of the total scattered 
power. 

d) Total differential cross sections. 

The cross section cr {6. ) per unit volume can be integrated analytically or 
even computed directly without considering the frequency domain (e.g. Hagfors 
1961, SaJpeter 1960, Kahn 1959, Renan 1960). The result for backscatter is: 




1 +a2 
1 + 2a2 


V.21 


This expression gives the result obtained by Gordon (1958) for a < < i ^nd 
reduces to 


a = cr 

e 



V.22 


for a > > 1 . 

It must be noticed that the power contained in the plasma line being pro- 
portional to 1/(2 a^), the expression V.21 gives in fact the power contained In 
the ion spectrum for a > > 1 , 
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e) Conclusion: 

It is worth noting that, fortunately, the problem of incoherent scatter was 
approached on the basis of random electron thermal fluctuations. Indeed the 
current theories of collective oscillations in the plasma at that time (e.g. Pines 
and Bohm 1952) would have led to the predictions of the plasma lines only and 
to a correspondingly negligible small scattered power, since the ions were then 
considered as a uniform background. 

The next sections concern the inclusion of magnetic field, unequal ion and 
electron tmnperatures, collisions, drifts of the electrons or ions and the effects 
of photoelectrons, 

VI. Case of Unequal Ion and Electron Temperatures 

The problem of nonthermal equilibrium is of considerable interest for the 
ionospheric plasma. Indeed cascading of energy from ttie photo-electrons to the 
thermal electrons, to the ions and finally to the neutrals implies that generally 
electrons, neutrals and ions have different temperatures. The extension of the 
computations of the previous section to the case of nonthermal equilibrium 
between the ions and the electrons was imdertaken by Salpeter <1960, 1963), 
Eejer 1961, Renau, Camnitz and Flood (1962), Rosenbluth and Rostoker (1962), 
Farley (1966). 

The results obtained by Renau et al. (1962) differed from those of the other 
authors. The disagreement was due to an improper extension of the Nyquist 
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theorem made by Renau, et al. (1962). Farley (1966) later indicated the valid 
extension of the Nyquist theorem to the case of nontherznal equilibrium between 
the ions and the electrons. The approach of Farley is ^ain outlined here. 

For a system in quasi equilibrium (Farley 1970), in viiich tiie electron and 
ion t^nperatures are not equal the equivalent electrical circuit is one containing 
impedances at different temperatures; the noise output is such that: 


<lE (<^)j ^ dw = 1 2^ [Z„ (co)] VI.l 

n 

where Zn is the impedance, and E the open circuit voltage. The corres- 
ponding expression applying to the incoherent scatter problem is (Farley 1966); 





K, [TA * T.Z.] 


VI.2 


vdiere 


and 

Z - 

e *ekk 

As a consequence the generalized expression for the power spectrum of the 
electron density fluctuations is; 
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where 

" KT^.i^o k2 

in the following a is defined as a = . While the electron spectrum is un- 

changed, the shape of the ion spectrum undergoes large changes when Te/Ti 
varies as illustrated in Figure 4 for a >> 1. For negligibly small values of 
the ratio Te/Ti, the ion spectrum tends to the gaussian velocity distribution of 
the ions. Indeed for the range of velocities of interest the momentum equation 
for the electrons reduces to a balance between the electron pressure gradient 
gnri the Space charge electric field. The pressure becoming negligible the space 
charge electric field tends to zero and the ion fluctuations are those of non- 
interactir^ particles. 

Increasing the ratio of Te/Ti has two effects: the first is to increase the 
space charge electric field thus reducing the ion fluctuations, this corresponds 
to the general decrease of the amplitude of the power spectrum for increasing 
Te/Ti. 
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The second effect is to enhance the pseudo acoustic waves, indeed the phase 
velocity of the waves increases as (Te+ Ti) and tiierefore less and less ions 
interact strongly with the waves when Te increases, thus reducing the damping 
of such waves. However further increases In the electron temperature lead to 
a damping of the pseudo sound waves dominated by the electrons rather than the 
ions (Farley 1966) with the net effect of increasing the total scattered power. 

For moderate values of Te/Ti the total differential cross section has been 
computed to be (Buneman (1962), Farley (1966): 


cr = cr n-, 

e 0 


a" 


1 


VIA 


1 + a2 (a2 + (T/Ti) + 1) (1 + a^) 


which for large values of a, yields: 


1 + T/T. 


VI.5 


This ejqjression would be valid if the ratio of the ion to the electron masses 
was infinite. 

The actual variation of o- with T^/T. was computed by Moorecroft (1963). '' 

+ -f 

Figure 5 shows such a variation for 0 and H , for n^ equal to unity and for 
a > > 1. Also shown on the figure is the variation given by the approximate 
expression VI.5 which is a limit for an infinite ratio of ion to electron masses. 

The transition between a spectrum dominated by random ion thermal fluc- 
tuations (T^ /r. < < 1) to a spectrum of weakly damped pseudo sound waves 
corresponds to a minimum of the total scattered power. 
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The variation of the cross section with respect to l/a for various values 
of Tg /Ti and for n^ equal to unity is given on Figure 6 (After Moorcroft, 1963). 

The variation given by the approximate formula is also shown for Te/X - co . 
A monotonic increase of the cross section as a function of l/a is observed for 
y/p =1, corresponding to a transition from an ionic to an electronic spectrum. 
A large departure from the approximate formula is observed for T^ /T. i 15 
and for small values of l/a ; Indeed the approximate formula is equivalent to 
neglecting the damping of the pseudo sound waves by the electrons; for large 
values of T^ /T. the waves are in thermal equilibrium with the electrons rather 
with the ions because the former participate most in the damping (Farley 
1966). 

VII. Plasma containing several kinds of ions 

The case of a mixture of ions was studied by Buneman (1961) and Moorcroft 
(1964). It is easily shown that W(6». ) in eqn. VI.2 must be replaced by 

n 

•where P is the relative concentration of the n'^^' ion 



The effect of a transition from 0 to He is illustrated on Figure 7 (after 
Moorcroft 1964). 
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More generally for several types of ions with multiple charges W{d.^) must 
be replaced by 

n 

■w4iere Z „ is the charge number such that 

= 1 

n 

Vin. Effect of the magnetic field 

The effect of a uniform magnetic field on the incoherent scatter spectrum 
was studied by Laaspere (I960), Fejer (1961), Hagfors (1961), Farley, Dougherty 
and Barron (1961), Salpeter (1961a), Renan et al. (1961) Farley (1966). 

The calculation, once again, can be conducted through the use of the Nyquist 
theorem approach (Farley et al., 1961). The problem reveals significant com- 
plications because of off diagonal terms in Yi and Ye. However, it can be shown 
that except for values of <p (the angle between the magnetic field and the wave 
vector k) extremely close to tt/ 2 only the diagonal elements are important. 

This is equivalent to letting c - and is called the "longitadiual approximation" 
(Bernstein 1958, Farley et al. (1961). 

Therefore equation VI.2 is still valid if W(0) is replaced by 


w( 0 , 4>. = 1 - ie 


sin^ <i> sin^ l/2$t- 1/4 COS dt vni.i 
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■vrtiere is the normalized gyro frequency: 


mk \2KT; 


vin.2 


The integral is called the Gordeyev integral. 

It can be shown (Farley et al. 1961) that for the ionospheric case no change 
in the spectrum is to be expected from the presence of the magnetic field for 
d, < 85“. For d > 85° the dispersion relation has new roots corresponding 
approximately to integral multiples of the ion or electron gyro frequencies 
(0 « n0» * e ^ n$ ) and called Bernstein modes (Bernstein, 1958). They 

correspond to particles which are periodically in the same phase relationship 
with the plasma waves. 

However due to the effect of the mass factor upon the gyrofrequency, the 
ion spectrum for practical purposes (a > > i) will only be affected by ion 
gyroresonances as illustrated on Figure 8 for d — 88°. The general shape of 
the spectrum is similar to the one obtained for ^ = 0° (or no magnetic field) 
with in addition a modulation at approximately the ion gyro frequency as originally 
suggested by Bowles (1959), When d 90° the gyroresonances peaks are more 
and more pronounced and eventually the spectrum becomes a spectrum of S. 
functions; however simultaneously all the power tends to be concentrated at the 
central frequency. 

It will be shown later that when ion collisions are included the modulation 
tends to be damped. 
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The total power scattered can be shown to be independent of the magnetic 
field for . The effect of the magnetic field for T^ T. can be shown 

to be equivalent to a change in the ratio of the mass of the electrons to the ions. 
The only quantity of interest is 



This is due to the fact that the electrons are stroi^ly bound to the magnetic 
field, a > > 1 corresponds to wavelengths v^ich are large both with respect to 
the Debye length and to the electron Larmor radius; as a result only the field 
aligned component of the electron velocity plays a part in the collective plasma 
oscillations: k • v^ must be replaced by k * v^^^ where Vgn is the field 
aligned electron velocity. This is equivalent to changing m]^^ into m^^ sec 0. 
(This change of mass obviously only applies to the determination of the plasma 
oscillations not to the scattering of electromagnetic waves). 

Figure 9 shows the variation of the normalized cross section as a function 
of T^ /T. for various values of y (after Farley (1966), 

The interpretation given by Farley is that sound waves get into thermal 
equilibrium with the particle population which contributes most to their Landau 
damping. The amplitude of the wave is then proportional to T„/ (T^ + T. ) in 
which T„ is the temperature of the population which contributes most to the 
damping. For small values of T^ /Tj , is equal to T. and for large values 
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of T /T T is equal to T . This explains why for both large or small values 
of Tg/Ti the ratio T„/T^ + T^ tends to 1. The transition corresponds to a 
minimum and depends upon the ratio of the effective masses: the larger the 
effective electron mass, the smaller the values of T^ /T^ for which the Landau 
of the sound waves is dominated by the electrons. 

The location of the plasma line is also modified by a uniform magnetic 
field (Salpeter, 1961b; Perkins and Salpeter 1965): 

0 ^= 0 ^ + ^ sin^ 4 > vni.3 

provided O-l << where is the electron gyrofrequency. In normalized 
units eqn. VIII.3 can be written as; 

6^2 vni.4 

^ er 2 ® 

Salpeter (1961b) has shown that an additional root of the dispersion relation is 
associated with the presence of the magnetic field (see also Perkins (1967): 

col - S vni.5 

in the ionosphere o> is generally larger than so that for <p - tt/2 

a,, ^ 

which is the geometric mean gyrofrequency. An additional line therefore appears 
in the incoherent scatter spectrum. 
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The effect of the magnetic field has been considered in this section for the 
case of practical interest for the ionosphere. However a number of other situ- 
ations might be considered according to the relative amplitudes of and the 
ion and electron Larmor radii. (See Farley et al. 1961) for example). 

IX. Effect of bulk motions of the plasma or of one type of particles. 

a. Motion of the whole plasma = v . 0 

It is straightforward to show that the v^ole ion spectrum is not changed in 
shape but is simply doppler shifted by So; = - k ‘ 

b. relative motion of the electrons with respect to the ions. 

This case has been studied by Rosenbluth and Rostoker (1962) and Tamh (1962). 

The spectrum becomes asymmetrical: the peak corresponding to the pseudo 
sound waves propagating in the direction of the electron drift is enhanced, while 
the other peak is damped. Indeed the electrons tend to feed the pseudo sound 
waves propagating in their direction and to damp the ones propagating in the 
opposite direction. This effect, however, only takes place when tiie electron 
drift (Vj) is a significant fraction of the electron thermal velocity (v^^) as shown 
in Figure 10 for ions, T^ /T. = 2, and backscatter (after Lamb 1962). 

The peak grows rapidly when vj approaches corresponding to the two 
streams instability. 

X. Effect of Collisions 

The collisionless case corresponds to many situations encountered in the 
ionosphere. It ceased to be valid, however, everytime a collision frequency 
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becomes of the same order of magnitude as the width of some line in the spectrum: 
this is the case in the E region when the ion neutral collision frequency is not 
negligible as compared to the width of the ion spectrum. This is also the case 
when gyroresonances tend to develop with widths of an order of magnitude close 
to the ion-ion collision frequency. And also, it is the case when the Landau 
damping of the plasma line becomes very small, electron-ion collisions must 
then be taken into account. 

The effect of collisions has been studied by Dougherty (1963), Cohen (1963), 
Dougherty and Farley (1963), Taylor and Comisar (1963), Waldteufel (1963, 1965), 
Tanenbaum (1968), Seasholtz and Tanenbaum (1969), Hagfors (1970), Hagfors and 
Brockeiman (1971). 

a. Ion-neutral collisions. 

The complexity of the problem is increased considerably when collisions 
are included. Indeed the simple relaxation term generally used for the Boltzmann 
equation is not convenient here since the scattering properties depend sensitively 
upon the detailed kinetic behavior of the plasma (Hagfors 1970). 

Efforts have been made to include phenomenologically collision terms 
which would conserve some of the physical invariants such as particle number 
density, momentum, energy, thus following the procedure of Bhatnagar-Gtross- 
Krook (1954) who proposed a term preserving charged particles through collisions. 

Figure 11 (after Dougherty and Farley 1963) shows the effect of collisions 
on 'he ion spectinim using a term preservii^ charged particles through collisions. 
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\p is a. normalized collision frequency; 



For increasing ion neutral collision frequencies the ion spectrum tends to 
become a Lorentz curve of width proportional to 

Hagfors (1970) has described a new approach to the problem, in which 
random forces exerted by the collisions with the neutrals are included among 
the forces which determine the unperturbed motion of the charged particles. 
The problem is then "to determine a transition probability which describes the 
diffusion of a charged particle through position velocity space" (Hagfors 1970). 
Hagfors and Brockelman (1971) applied it to hard spheres elastic collisions. 
The gross features are identical to the one of Figure 10, however the central 
depression of the spectrum does not disappear as rapidly when */'i increases. 

b. Ion-Ion collisions. 

In the F region of the ionosphere the frequency of ion-ion 90® collisions is 
generally smaller than the Ion gyro frequency. However small angle deviations 
through collisions are sufficient to affect significantly the ion trajectories over 
one gjroperiod. As a consequence gyroresonances are completely destroyed 
over most of the F region, 

c. Electron- Ion collisions 

It was mentioned earlier that Landau damping tends to bring plasma waves 
in thermal equilibrium with the plasma. This holds for the plasma lines. 
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however when the Landau damping becomes small because of the small number 
of electrons participating in it, damping throu^ electron-ion collisions becomes 
dominant. In these circumstances this last mechanism is responsible for bring- 
ing the plasma waves in thermal equilibrium with the plasma (Perkins and 
Salpeter (1965)). 

The total intensity of the plasma lines is consequently unchanged but their 
width is now determined by collisional damping. 

Interesting situations develop when the excitation of the plasma line by 
photoelectrons is taken into account. This is ttie subject of the next section. 

XI. Plasma lines enhanced by the photoelectrons 

Plasma lines enhancement by photoelectrons have been studied by Perkins 
and Salpeter (1965), Perkins et al. (1965), Yngvesson and Perkins (1968). 

When a tail of photoelectrons is added to the thermal electron Maxwellian 
distribution, the plasma lines reach thermal equilibrium with either the theimal 
plasma or the photoelectron population; in the latter case the equivalent tem- 
perature is much higher and the intensity of the line is correspondingly much 
higher. The particle population which determines the intensity of the wave 
depends upon the relative strengths of the Landau dampings by the thermal 
electrons or the photoelectrons and of the electron ion collisions damping; if 
either the thermal electron Landau damping or the electron-ion collision damp- 
ing dominates, the plasma line is in thermal equilibrium with the thermal plasma; 
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while, if the Landau damping by photoelectrons dominates, the plasma line is in 


thermal equilibrium with the photoelectron population. 

Perkins and Salpeter (1965) derived the following formula for the enhance- 
ment of the plasma line intensity in a plasma without magnetic field with respect 
to its level when photoelectrons are absent: 

e 

where f is the one dimensional velocity distribution of the thermal electrons 
in the k direction, fp is the photoelectron velocity distribution in the same 
direction, x describes the electron ion collision damping and is the wave 
phase velocity. 




±1 

6tt^ 


(2^) 


Log^(47rnoX3) 


XI.2 


The dominant term in eqn. XI.l is the thermal electron damping for low 
values of , the collisional damping for large values of and the photoelectron 
damping for intermediate values of v^. 

Therefore the plasma line is enhanced only over a limited range of values 
of or in other words over a limited range of resonant frequencies. 

The effect of the magnetic field is to considerably increase the Landau 
damping by the thermal electrons and therefore to generally decrease the range 
of frequencies over which the plasma line is enhanced. 
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Salpeter (1961b) has shown that the effect of the magnetic field can be taken 
into account by modifying the velocity distribution of the thermal electrons in 


the following way; 

, 1/2 
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XI.3 


where I is the Bessel function of imaginary argument. When 0 varies from 0° 

TO 

towards 90“, the frequency range dominated by thermal electron Landau damping 
increases steadily. However, when 90“ is approached the Landau damping by 
thermal electrons becomes strongly modulated with maxima at frequencies equal 
to multiples of the electron gyrofrequency. Eventually the Landau damping by 
thermal electrons disappears when 0 reaches 90°. 

These features are illustrated on Figure 12 which shows the frequency 
ranges of plasma line enhancements for various angles 0 and wavelengths after 
Yngvesson and Perkins (1968) and for a given f^ distribution. The enhancement 
is expressed in terms of plasma temperatures, which implies that when fully 
enhanced the line is in thermal equilibrium with the local photoelectron velocity 
distribution characterized by an apparent temperature T^ . Changing the wave- 
length has the effect of changing the phase velocity associated with a given 
plasma frequency; as a consequence the range of frequencies over which the 
■plasma line is enhanced depends upon the wavelength. In turn, since for practical 
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purposes the plasma frequency in the ionosphere covers a limited range the 
overlap of the two frequency ranges will be more or less favorable according 
to the wavelei^th used. 

For a wavelength X of .70 m there is a good overlap and one can see the 
three regions: the region dominated by thermal electron Landau damping for 
low frequencies and its eiqiansion for increasing values of the region where 
the line is enhanced because of the predominance of the Landau damping by the 
photo electrons for intermediate values of the plasma frequency, and finally the 
region dominated by electron-ion collisional damping. For = .23 m, the 
region dominated by thermal electron damping has expanded, while for ^ = 200 m 
most of the range corresponds to collisional damping. The modulation at the 
electron gyrofrequency is quite pronounced in all cases for = 80". 

Xn. Conclusion 

The simple concept of independent free electrons initiated the theory of 
incoherent scatter. The problem became much more complex \^en it was 
realized that collective interactions had to be taken Into account for various 
situations. 

The nature of the complexity being associated with some of the essential 
properties of the plasma, the elaboration of the theory became more and more 
rewarding for both the theorist and the experimenter. 

In particular in practically all instances, every effort made by the theorists 
to include a new parameter in the theory was accompanied by an immediate 
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experimental verification, the experimental work thus helping considerably the 
theoretical work. 

In turn the addition of new parameters was to permit the e>perim enters 
a more and more complete diagnosis of the state of the ionosphere. The 
simultaneous determination of 4 to 5 ionospheric parameters has certainly 
made incoherent scatter the most powerful ground based mean for studying the 
dynamics and thermodynamics of the ionosphere and the thermosphere: a 
considerable contribution to these tields has already made by incoherent scatter 
soimding (e.g. reviews by Evans (1969, 1964) and by Farley (1970b), Incoherent 
scatter sounding using laser beams is also more and more used for laboratory 
plasma diagnosis. 

The future of incoherent scatter sounding of the ionosphere involves, In 
addition to observations from the existing facilities, important projects in 
North America and Northern Europe and even a possible implementation of an 
incoherent scatter radar on board the Space Shuttle Space Laboratory. In all 
instances the future of incoherent scatter for the study of the ionosphere is to 
result from a joint effort of scientists from many countries. 
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GEOMETRY OF INCOHERENT SCATTERING 
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TRANSM ITTER RECEIVER 

Figure 1. Geometry of incoherent scattering. 
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NORMALIZED DOPPLER SHIFT 9 , 


Figure 2. Normalized spectrum of the ionic component for a » 1 and T = T. 
(full line). Ion velocity distribution (dashed line). The normalized doppler 
shift is defined by 




= V 
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(after Farley, 1970a). 
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POWER cr(ee)v^/nor 


COMPLETE SPECTRA FOR VARIOUS VALUES OF a 



Fi^e 3. Normalized complete spectra for different values of a = , 

(Tg = T. ). The normalized doppler shift is 

.3 r m I ^''2 

e 

® k 2KT 

(after Hagfors, 1961). 
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IONIC COMPONENT VS. T./Tj 
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Figure 4. Normalized half spectra of the ionic component for a >> i and for 
different values of Tg/Ti. {After Fejer, 1961), 
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Figure 5. Variation of the ratio of the total differential cross section for a » 1 to the total differ- 
ential cross section for independent electrons as a function to /T. (after Moorcrost, 1963)* The 
dashed curve corresponds to the approximate expression VI.5 given by Buneman (1962). 
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Figure 6. Variation of the ratio of the total differential cross section <r to 
the total differential cross section for independent electrons as a function 
of a, for H+ and for different values of T /T. (after Moorcroft, 1963). 

The dashed line corresponds to the approximate expression given by 
Buneman (1962), 
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RELATIVE POWER 



Figure 7. Half spectra of the ionic component (normalized to a value of 1 for zero doppler 
shift) for Tg/Tj =1.5, a » 1 and for different mixtures of 0+ and He+, N2/N is the relative 
concentration He+. (After Moor croft, 1964). 
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ION GYRO-RESONANCE LINES 



Figure 8. Normalized half spectra of the ionic component for 
Tg = Tj , a 0+ ions, A. ~ R, (the Lamor radius of the 
ions), and for k either parallel to the magnetic field = 0% or 
no magnetic field) or nearly perpendicular to it (<?b = 88°). Ion- 
Ion collisions have been assumed to be negligible. (After Farley 
et al., 1961). 
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EFFECT OF 

ELECTRON DRIFT ON IONIC COMPONENT 

Te/Ti = 2.0 (3=102 



Figure 10. Normalized spectra of the ionic component for 
T /T. ~ 2.0 a =10^, H+ ions and different electron drifts 
antiparallel to k. v^/v^ is the ratio of the drift velocity to 
the electron thermal velocity. (After Lamb, 1962). 
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NORMALIZED DOPPLER SHIFT 0j 


Figure 11. Normalized half spectra of the ionic component for 
T =T. ,a>>l and for different normalized ion neutral collisions 

e X * 

frequencies 



(After Dougherty and Farley, 1963). 
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PLASMA FREQUENCY (Mhz) 


Figure 12. Plasma line enhancements for various wave-lengths 
and angles 4 > between the k vector and the ms^etic field as func- 
tions of the plasma frequency. The photoelectron velocity dis- 
tribution fp was deduced from measurements and was chosen 
equal to 2000 K (after Yngvesson and Perldns, 1968), 
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